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We study the real-time dynamics of the order parameter (cr(i)) in the Ising field theory after 
a quench in the fermion mass, which corresponds to a quench in the transverse field of the cor- 
responding transverse field Ising chain. We focus on quenches within the ordered phase. The 
long-time behaviour is obtained analytically by a resummation of the leading divergent terms in 
■".."_ a form-factor expansion for {o(t)). Our main result is the development of a method for treating 

, ■ divergences associated with working directly in the field theory limit. We recover the scaling limit 

^— ^ ■ of the corresponding result by Calabrese et al. [Phys. Rev. Lett. 106, 227203 (2011)], which was 

-^^ ' obtained for the lattice model. Our formalism generalizes to integrable quantum quenches in other 

, integrable models. 

"*^ ! I. INTRODUCTION 

o; 

^^ ■ Recent experimental advances have made it possible to study the non-equilibrium dynamics of trapped cold atomic 

' gases^i^. A key feature of these systems is that they are only weakly coupled to their environments, which makes 

it possible to study non-equilibrium dynamics in essentially isolated systems. This has led to an intense theoretical 

effort to address fundamental questions^ regarding the non-equilibrium dynamics of many-body quantum systems. 

One issue of particular interest concerns the role played by dimensionality and conservation laws. As shown by the 

^ ■ "quantum Newton's cradle" experiments of Kinoshita et al.^, quasi-one dimensional condensates exhibit behaviour 

that is dramatically different from two and three dimensional ones. In particular, it was observed that the late-time 

behaviour cannot be described in terms of an effective temperature: the systems does not "thermalize'"^. One possible 

Ch ' explanation^ for this behaviour is that the experimental system is close to being integrable. This has engendered a 

1-^ . vigorous research effort aimed at clarifying the role played by quantum integrability in the late time and stationary 

(— I ' state behaviour in non-equilibrium dynamics^^— . 

O . A simple and attractive way of inducing non-equilibrium evolution is by means of a quantum quench. One prepares 

O ■ a system in the ground state of a given Hamiltonian H{h()), where Hq is an experimentally tuneable parameter such 

[ as a magnetic field or an interaction strength. At time t = the parameter ho is then changed instantaneously 

p^ ■ from ho to h, and at subsequent times the system evolves according to the quantum dynamics induced by the new 

^ ■ Hamiltonian H{h). One of the main models studied in the context of quantum quenches has been the transverse field 

^D . Ising chain^*^"^^. This is on the one hand because the model has a free fermion representation which makes analytical 

2~ ■ progress possible. On the other hand the model is the simplest paradigm of a quantum phase transition and therefore 

is an ideal testing ground for questions relating to non-equilibrium evolution in the vicinity of quantum critical points. 

The stationary and late-time behaviour of correlation functions in the transverse field Ising chain after a quantum 

CO I quench has recently been determined analytically by Calabrese, Essler and Fagotti^°. 

^— ^ ■ Analyzing quantum quenches in interacting integrable models is difficulti'^ and remains a largely open challenge, 

although important progress has been made by combining numerical and integrable model techniques*. A special role 

is played by integrable quenches in integrable quantum field theories. These are characterized as follows. As shown 

^ ■ by Calabrese and Cardy— the quench problem can be mapped to an equivalent theory defined in a strip geometry. 

k>( . The initial state plays the role of a boundary condition, and for an integrable quench this boundary condition does 

j_j ■ not spoil the integrability of the theory. Hence, for these special initial states one can use methods of integrable 

C^ ] quantum field theor y ^"^i^^ with boundaries'^^ to analyze the time evolution of observables. An important step in this 

direction was taken by Fioretto and Mussardo^, who considered the stationary state behaviour of one-point functions 

in integrable quenches and in particular in the Ising field theory. A serious complication that arises in the field theory 

limit is that singularities associated with kinematical poles appear. This problem is particularly acute for two-point 

correlators and is closely related to the one encountered when calculating finite-temperature dynamical correlation 

functions in integrable modelaiSr— . To date two general ways of dealing with these singularities have been developed. 

The firsiii^r— is to use a finite- volume regularization for matrix elements^^, while the second is a subtraction scheme 

that works directly in the infinite volume^^. The aim of the present work is to apply these methods to the problem 

of quench dynamics in the ordered phase of the Ising field theory. Regulating the theory in a finite volume reduces 
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all calculations to a particular limit of the analysis for the lattice Ising modeliSi^l and we therefore do not report any 
details here. We focus on the infinite-volume regularization proposed in Ref. [17[ and apply it to the quench problem 
at late, finite times. This requires a significant generalization of the regularization procedure, which constitutes the 
main result reported here. In forthcoming work we will apply this method to a quench in the sine-Gordon model. 

The outline of this paper is as follows: In Section |ll] we introduce the Ising field theory as the scaling limit of the 
transverse field Ising chain. In Sec. lIIII we discuss quenches in the fermion mass of the field theory, which corresponds 
to a quench in the transverse field of the related Ising chain. In Sec. lIVI we develop a new method to calculate the time 
evolution of correlation functions in integrable field theories, which constitutes the main result of our work. In Sec. |V] 
we apply this method to the one-point function of the order parameter field, which relaxes exponentially to zero as 
shown in Eqs. (H51) and (H71) . In Sec. |VT we discuss the relation of quenches in the fermion mass to the extrapolation 
time regularization introduced in Ref. y| and conclude in Sec. I VIII Technical details of the derivations have been 
moved to the appendices. 

II. QUANTUM ISING CHAIN 

We start with the transverse field Ising chain 

H,,,, = -jJ2{a!a!+, + h<jf). (1) 

i 

Here a^ and a^ are the Pauli matrices and J is the exchange energy. The Hamiltonian ([T]) is invariant under the 
Z2-transformation af — >■ erf, erf — )■ —erf. For h < I this symmetry is broken, the order parameter field erf takes a 
non-zero expectation value, and the ground state is two-fold degenerate. On the other hand, for h > 1 the system 
possesses a unique ground state and the expectation value of the order parameter field vanishes. The two phases are 
separated by a quantum critical point at h = I. At small deviations from criticality, |/i — 1| <i; 1, one can pass to the 
scaling limit2^ {uq is the lattice spacing) 

J ^ oo, h^l, Co -> 0, (2) 

while keeping fixed both the gap M and the velocity v 

2J\l-h\^M, 2Jao = v. (3) 

The order parameter in the scaling limit must be defined as 

aix)^il-h^r^cTl (4) 

where x = nao- It is customary to choose the normalization of the field (j{x) such that 

hm {0\a{x)a{0)\0} = ^, (5) 

which implies 

^z ^ 2l/24el/8^-3/2„l/8^(^)^ (g) 

with Glaisher's constant 

yl= 1.28242712910062... (7) 

The Hamiltonian in the scaling limit reads 

dx 



^ ■ 2. 
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where ip and V' are the two components of a Majorana fermion. The model (|H1) is conformally invariant at the critical 
point M = (see for example Ref. [23]). In the ordered phase, which we will consider throughout this paper, the 
mass is positive. 

An important notion is the mutual semi-locality of operator o^^'"'^^'^'^'^^ . This is most easily established by defin- 
ing complex coordinates z = T + ix, z = t — ix and then considering the operator product Oi(r, x) ©2(0, 0) = 



Oi(z, z) 02(0,0). If we take Oi counterclockwise around O2 in the plane, i.e. we perform the analytic continuation 
z — ^ e'^'^^z, z — >■ e~^'^'z, the operators Oi and O2 are said to be mutually semi-local if 

Oiie^^'z, e-2-z) 02(0, 0) = /oi02 Oi(z, z) 02(0, 0). (9) 

The phase I01O2 is called the semi-locality factor. The two fields are mutually local if I01O2 = 1- Semi-locality is the 
mildest form of non-locality, in general the right-hand side of © may be more complicated. The mutual semi-locality 
factor of the spin and disorder operators can be extracted from their operator product expansio n^^i^^ 

a(z, z) ^,{0, 0) ^ ^1 [e'-/4V^^(0) + e-'^^'V^m] ■ (10) 

This implies that when taking a once around /i one obtains an extra minus sign, i.e. l^^ = — 1. In the same way one 
finds l^fj_ = Ijji^ = l^a = l^a — ~^- On the other hand, the disorder field /i is local with respect to itself. 

We use the disorder field /x as fundamental field creating the excitations. This implies that the fundamental 
excitations are viewed as bosons. We denote the corresponding annihilation and creation operators by A{6) and A^{0) 
respectively. They fulfil the Faddeev-Zamolodchikov algebra3^ 

A{eM{02) = sA[e2)A{ei), 
A\e^)A\e2) = SA\e2)A\6{), (ii) 

A{ei)A\e2) = 2TT5{6i-62) + SA\e2)A{ei), 

with the scattering matrix S" = — 1. The basis of scattering states can now be constructed by 

|0i,...,0„)=At(0i)...At(0„)|O), (12) 

where the vacuum state |0) is defined by A{d) |0) = 0. The energy and momentum of the scattering states are 
expressed in terms of the rapidities 9i as 

£; = M^cosh6'i, P=— ^sinh6l,. (13) 

In this article we study the one-point function of the order parameter field a. The relevant matrix elements (form 
factors) in the ordered phase are given by^iiS^iSl 

2" B - 

/(0i,...,02„)-(O|a|0i,...,^2n)=i"fT n tanh^— ^, (14) 






where 



^ 2l/12g-l/8_43/2 I'M^ . (15) 



III. QUENCH IN THE FERMION MASS 

We now consider a sudden change of the transverse field in ^ at time t = from ho to h. This quench has 
been studied previously by several authors^iiiii^. Most importantly, in Refs. [10||21 | the time evolution of both the 



one-point and two-point function of the order parameter after a quench was determined analytically. One of the 
methods developed in Refs. |10ll2l] | is a form-factor approach for the lattice model. In the following we consider the 
time evolution of the one-point function directly in the scaling limit ([2]), ([3]). As we have mentioned before, our key 
objective is to generalize the form-factor approach to quantum field theories in order to analyze integrable quenches 
in interacting systems such as the sine-Gordon model. However, a second interesting issue is related to commutativity 
of limits: a priori it is unknown whether a quench in the scaling limit is the same as the scaling limit of a quench. 
We will come back to this question in Sees. IVFI and lVIII In the following we resolve this question for the particular 
case of the one-point function of the order parameter in the ferromagnetic phase of the Ising model. 



In the field theory ([5]) the quench in the transverse field corresponds to a quench in the fermion mass, i.e. at time 
t = we switch from Mq to M. The time evolution for t > is governed by ([5]), while the initial state can be 
expressed in terms of the eigenstates of ([8]) aai^iH 

\^o) - cxp ^j^ ^Xq(0At(-0^t(0) |0) , (16) 



where 



^q(C) = itan 



— arctan(sinh^) — — arctan — -sinh^ 



= i^q(e). (17) 



We note that the quench matrix satisfies -ft'q(^) = S K^{—^) = — ii'q(— ^) and that Kq{,0 G R for ^ G M. Furthermore, 
for any finite initial mass, Mq < oo, the integral 
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(18) 



is convergent. We note the similarity of the initial state (ITCl) with the boundary statc^^ introduced in the context of 
integrable field theories with boundaries, which can be used to study the physical properties of systems with defects 
or impuritie o^^i^^ . Starting from the initial state (fTB|) we calculate the time evolution of the one-point function of the 
order parameter field, 

IV. METHOD 

The strategy to calculate the one-point function p^ after the quench is as follows: (i) We formally expand^ i^^i-'^^ 
the numerator and denominator in powers of the quench matrix Xq. (ii) We evaluate each term in these expansions 
using a combined approach based on a regularization of the appearing form factors following Smirnov'^^^ as well as 
the K-regularization recently introduced in the study of dynamical correlation functions at finite temperatures^^, (iii) 
In the resulting expression the singularities in the numerator and denominator, which are due to the infinite volume 
of the model (H]), cancel each other. In particular we show by explicit calculation up to ©(i^T*) that this procedure 
yields well defined results which agree with a finite-volume regularization. (iv) Finally the resulting series has to be 
resummed"^? in order to obtain a well-defined long-time limit. The calculation of two-point functions follows the same 
lines, although the explicit expressions become considerably more complicated. 

A. Formal expansion 

The first step in the calculation of (jl9|) is the formal expansion of both the numerator and the denominator in 
powers of the quench matrix A'q. This expansion yields for the numerator 

X (ei, -el, . . . , e™, -Ck |-?«, Cn, • ■ • , -Ci, Ci) (20) 

OO 

= J2 C^rnMt)- (21) 

m.n—Q 

Note that the indices 2m and 2n correspond to the number of particles originating from the left and right initial state 
respectively. Similarly, the expansion of the denominator reads 

(*o|*o) - 1 + 1:/°° ^S(2.?»" t!(2.?»" n K,{Q*K,i^.) {e^,-e,,.-., C. -O - C», C», ■ . ■ , -Ci, 6) (22) 

OO 

- 51 ^2„. (23) 
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The normalization in (J19p can therefore be formally expanded in the following way 

We note that ([21]) is merely used for defining linked clusters, i.e. identifying the parts of the numerator in ([T5]) that 
diverge in the infinite volume. 

B. Regularization procedure 

The matrix elements in the terms C2m,2n{t) and Z2n possess kinematical poles whenever ^^ = ^j and therefore have 
to be regularized. Following SmirnovM we proceed as follows: Let A denote a set of one-particle excitations and Ai 
and A2 a partition of A. The scattering matrix arising from the commutations necessary to rewrite \A) as IA2A1) is 
denoted by SaAi, i-e. \A) — SaAi \A2A1) = Saa2 |^i^2)- If A and B denote two sets of one-particle excitations, the 
form factors of a readi^i^ 

{A\a\B) ^ J2 diB2)SAA,SB,B {A2\B2) {Ai+iO\a\Bi), (25) 

A=AiUA2 
B=BiUB2 

where the sum is over all possible ways to break the sets A — A1UA2 and B — B1UB2 into subsets. The scalar products 
(^2! B2) as well as the corresponding terms in the Z2n& are easily evaluated using the Faddeev-Zamolodchikov algebra 
([TTI) . The factor d{A) is present by virtue of the semi- locality of the spin operator with respect to the fundamental 
field and is given by 

d{A) = (-1)"(-^), (26) 

where n{A) denotes the number of elements in A. Using (j25l) the poles in the form factors have been shifted away 
from the real axis. The terms (A2I B2) correspond to the disconnected pieces of the form factors. As all rapidities in 
the remaining matrix elements are distinct, they can be evaluated using the crossing relation 

{e[ + iO, . . . , C + iOk 1^1, . . . , 0„) = f{e[ + Itt + iryi, . . . , C +m + It^™, 0i, . . . , 0„), (27) 

where 7?^ — ?► 0+. We note that one can also shift the rapidities in the set Ai to the lower half plane, which results^^ 
in diflFerent scattering and phase factors in (P5|) but leaves the final result unchanged. 

It is clear that the right-hand side of (P5]) may still contain divergences due to the intertwining of particles with 
rapidities ^i and —^i in the initial state (J16l) . These divergences are a consequence of working in the infinite volume and 
have to be canceled against similar divergences originating from the norm of the initial state (1^51) . In order to exhibit 
these cancellations we need to identify these divergences explicitly. To this end we use the K-regularization scheme 
recently introduced in the study of finite-temperature correlation functionaii. For each pair of rapidities {— ^i, ^i} in 
the ket states we introduce an auxiliary real parameter Hi to shift the rapidities away from the singularities. The 
resulting expressions have to be understood as generalized functions of the auxiliary variables k^. In order to exhibit 
the cancellations of terms in the Lehmann representation of (J19p that diverge in the infinite volume we define a smooth 
function P{k) which is strongly peaked around k = and satisfies 

P{0)^L, jdKP{K)^l. (28) 

Here L can be thought of as the length of the system in the finite- volume regularization (see App. [FJ. One possible 
choice is P{k) — L e"'^^ " . Using this regularization scheme the first non-trivial term in the expansion ((23l) reads 



dKP{K)S{-2K) d^Kq{^ + K)*Kq{0 

^inax{0, — /t} 

■niax{0, — k} 



(29) 



//•max-^u, — K| 
dKP{K)5{-2n) I d^Kq{^^-K)*Kq{^) (30) 

/"GO 

/ d^\KM'. (31) 

Jo 

The infinite-volume divergence is now clearly exhibited and (I3ip facilitates comparison with the finite-volume reg- 
ularization (see App. IfJ. Further examples for the application of (1251) and the k- regularization are presented in 

App.EHl 
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C. Cancellation of singularities 

Using the K-regularization as described in the previous section all terms in the expansions (j2ip and (1231) are finite 
but may contain terms ex L*^, or equivalently ex (5(— 2k,;), that diverge in the infinite- volume limit. However, when we 
consider the one-point function (J19p given by their quotient and expand again in powers of Kq, 

Eoo ^ /,N oo 

{a(t)) = "'^-c^ „' '' ^ J2 D2rr..2nit), (32) 

^-^n—v rrt,n—0 

all terms ex L'^ with fc > 1 cancel each other and the remaining functions D2m,2n{t) are finite in the infinite- volume 
limit L — >■ oo. This can be thought of as a linked-cluster expansion and is analogous to the finite-temperature caseii. 

D. Resummation 

Performing the steps outlined in the proceeding sections we obtain the expansion (j32[) for which the infinite- volume 
limit can safely be performed. After taking this limit we can study the long-time behaviour of the one-point function. 
Doing so we observe that the leading contribution to the term D2m,2n(i) will grow as ex i" with the power a depending 
on the number of particles 2m and 2n. As we will show below, however, these divergences can be resummed leading 
to a well-defined long-time behaviour of the one-point function, which we present in Sec. IV Fl 

V. RESULTS 

In this section we present the results for the leading terms in the expansion (|32l) . We consider terms up to 0{K^) 
and concentrate on the dominant contributions in the long-time limit. Technical details of the derivation are presented 
in Adds. EHdI In Sec. IV Fl we present the final result for the long-time behaviour of the one-point function (J19l) after 
the resummation of the leading contributions in the series p2p . 

A. Terms in 0{K°) and 0{Kq) 

The terms up to linear order in Kq do not contain form factors possessing both incoming and outgoing particles. 
Thus there exist no kinematical poles, a regularization following the procedure discussed in Sec. lIVBJ is not necessary 
and we straightforwardly obtain 

-Doo = Coo = (0| a\0)^a (33) 

as well as 

D2oit) + Do2{t) = C2o{t) + Co2{t) =a r ^ i^q(0 tanhefe^*"*-^^? ^ p-2Afi*cosh5l (34) 

Here the real function Kq{^) was defined in (J17p . The long-time behaviour of this term is obtained by a stationary 
phase approximation, 

DMt)+Do2it)^^^[l--) \,,,^J\ Mt»l. (35) 

B. Terms in 0{K^) 

In this order there exist three terms. The first two originate from Cio{t) a-nd Co4(i), which do not possess kinematical 
poles. In the long-time limit we obtain D^oit) + Doi{t) ^ cos{AMt) / [Mt)^ which constitutes a sub-leading correction 
to ([35]). 

In contrast C22 contains a form factor possessing both incoming and outgoing particles and hence kinematical poles 
appear. Performing the calculation as outlined in Sees. IIVBI and I IV CI we obtain (see App. |B] for details of the 
derivation) 

D22{t) = C22{t)-Z2a^~aTt + D'^^{t), (36) 



r = — / rf^|/^q(OpsinhC, (37) 

Here the contour of integration 7_ lies in the lower half plane and can be explicitly parametrized by {0 < (j)o < 7''/4) 

7-(^)-(r"'^^ -^ ^-<-/°' (39) 

i (s - (po) - Kpo, 00 < s < oo. ^ ^ 

Clearly the first term in (l36l) dominates the long-time behaviour. It can be thought of as the second term in the 
expansion of tr e~^* in powers oi Kq, see Sec. lVFI The late-time behaviour of the second contribution psp is dominated 
by the region ^' « and ^ = —is with s ~ 0. Expanding the integrand and changing to polar coordinates then gives 

D',,it).-^^1^^^IM^, Mt»l. (40) 

We note that the behaviour of p6p is in agreement with the finite- volume regularization presented in App. |F1 

C. Terms in 0{K^) 

In this order the terms containing kinematical poles are C42(i) and C24(t). The calculation presented in App. \C\ 
yields 

Di2{t) - D24{tr - C42{t) - Z2 C20{t) = -T t D2o{t) + ..., (41) 

where the dots represent sub-leading terms that fall off at least as ~ l/(Mi) in the long-time limit. Again we find 
a term showing an explicit linear time dependence, which can be viewed as the second term in the expansion of 
D2o{t)e-^'. 

D. Terms in 0{K^) 

The calculation of D44{t) requires the introduction of two independent auxiliary parameters ki and K2 (see App. lDJl . 
The result is given by 

D^iit) = C44(t) - ^2^22(0 + [Zl - Zi)a = |(rt)2 -rtD'^^{t) + D'^^ + ..., (42) 



^44 - -a^t j^'^\Kq{^ 



'^^'-,(6)1^ 



§7^,(6)^ (coth^ ^ - tanh^ ^ - coth^ ^ + tanh^ ^) , (43) 

where the dots again represent sub-leading terms that fall off at least as '^ l/{Mt) in the long-time limit. The first 
term in (|42p can be viewed as the third term in the expansion of ae"'"*, while the second corresponds to the second 
term in the expansion of i?22(0 e '"*• On the other hand I?44 is independent of time and represents a correction of 
order K^ to Doo = ^- We further note that £'44(t) does not contain terms that are linear in t [as D'22{t) oc l/(Mi)]. 

E. Leading time dependence of higher-order terms 

Finally we argue in App. [E]that the leading term in the long-time behaviour of D2m,2m{t) and D2m+2.2m{t) are 
given by 

D2m,2m{t) = ^{-rtr + ..., (44) 

D2m+2,2m{t) = i?2™.2™+2 (t)* - ^(-^0'" D2o{t) + ..., (45) 

m] 
where the dots represent terms that grow at most as oc t™^^ for large times. 



F. Resummation and long-time behaviour 

As we have shown in the previous sections, the expansion p2p contains terms that grow at long times as powers of 
t. Hence, in order to obtain a weU-defined result in the long-time limit we have to resum these divergences. From the 
results presented in the previous sub-sections we deduce the leading long-time behaviour (Mt ^ 1) of the one-point 
function of a after a mass quench. 



{<j{t))=<j 



1 



a 

m 



1 - M/Mo cos{2Mt - tt/4) 



8v^ 



{Mt)y^ 



-,-rt 



(46) 



where the dots represent sub-leading corrections to the prefactor as well as terms that decay faster than e 
respectively. The relaxation rate is given by 



-rt 



2M f°° 

r(M, Mo) = / d^ |Xq(e)P sinh^ + 0(X, 

I" Jo 



(47) 



We have determined the relaxation rate up to order K^. The fact that there is no contribution in 0{K^) follows from 
the absence of terms linear in t in 1)44 (t). This finding is in complete agreement with the corresponding result for 
the lattice model^°. The consistent calculation of the relaxation rate in 0{K^) would require the derivation of the 
contributions to DeQ{t) that grow linearly in time. 

The 1 and the l/(Mi)^/^ term in the prefactor of the exponential in (|15|) have been established by considering 
particular contributions to all orders and showing that they exponentiate [see (j44l) and (|45)) 1. A detailed discussion 



of this point for the lattice Ising chain is given in Ref. (21|. On the other hand, the '-^ l/{Mt) contribution in the 



prefactor is a conjecture based on our results for the leading contributions (in the expansion in powers of the quench 
matrix) in the Lehmann representation of the one-point function. From (j40p and (j42l) we deduce 



(1 - M/MoY 
327r 



OiK^J 



(48) 



Finally we stress that our results agree with the scaling limit of a quench in the transverse field of the Ising chain in 
the ordered phas o^°i^^ . 



VI. EXTRAPOLATION TIME 



The AT-matrix for fixed boundary conditions in the Ising field theory is given byi^ 

A'fixod(C) 

We note that this is obtained as a limit of the quench AT-matrix (fT7|) 



i tanh — . 
2 



^fixcd(e) = , lim KgiO- 



(49) 



(50) 



In the quench problem a finite value of Mq is required to render rapidity integrals convergent at large energies. In 
particular, the decay rate F depends on Mq and diverges in the limit Mq -^ 00. For quenches in interacting integrable 
quantum field theories it is currently not known how to express a given initial state in terms of eigenstates of the post- 
quench Hamiltonianji^. An exception are initial states that correspond to integrable boundary conditions. In order 
to use this information in the context of quantum quenches, a prescription how to "regularize" the corresponding 
AT-matrices at large rapidities is required. Fioretto and Mussardo introduced an "extrapolation time" tq by the 
replacemenliS 



infixed (0 -^ itanh I e-2*^-« ™-h« ^ j^^^^f^^y 
Using this regularized A'-matrix to perform our calculations results in a decay rate 



2 A// r°° 

Tr, (M) = / d^ I AVo m' sinh f + 0(A:6 
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,)• 



(51) 



(52) 



By comparing (l52l) and (|47)) and requiring the decay rates to be equal Trg{M) — T{M, Mq) is it possible to relate the 
extrapolation time tq to the initial mass Mq. 



VII. CONCLUSIONS 

In this work we have considered the time evolution of the order parameter after a quantum quench of the mass 
in the Ising field theory. We have focussed on a quench within the ordered phase. We find exponential decay of 
the order parameter to zero (|46p . Our results agree with the scaling limit of a quantum quench performed in the 
ordered phase of the transverse field Ising chainiSiSi. Our main achievement is of technical nature: we have shown 
how to carry out calculations in the field theory limit. Here, unlike for the lattice model, additional divergences 
occur that need to be regulated appropriately. We have shown how to use techniques developed recently in the study 
of integrable quantum field theories at finite temperatures to overcome this problem. Our method generalizes to 
interacting integrable quantum field theories such as the sine-Gordon and 0{N) non-linear sigma models. This opens 
the door for analyzing quantum quenches in these theories, at least for particular classes of initial states related to 
integrable boundary conditions ("integrable quenches")®. Work on the sine-Gordon model is under way. 
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Appendix A: Calculation of {^o\ *o) 

In this appendix we evaluate the leading terms in the expansion 1^^ of the norm of the initial state |^o)- Obviously 
one has Zq = 1, while Z2 was already calculated in Sec. lIVBI with the result 



/"OO />oo 

/ dC\K^{0\'^S{-2K) d^\K^{0\'- (Al) 

Jo Jo 



In the last step we have reintroduced the auxiliary variable k. We will explicitly retain the auxiliary variables Ki 
throughout the appendices, but keep in mind that all expressions have to be understood as generalized functions of 
Ki as discussed in Sec. lIVBl As all expressions are multiplied by the strongly peaked functions P{Ki) we can drop all 
terms oc k" with n > 1. In contrast, all irregular terms ex S{Ki) as well as all divergent terms oc k" with n < —1 have 
to cancel when considering the expansion ([5^ . It is the purpose of these appendices to show by explicit evaluation up 
to 0{K^) that these terms indeed cancel each other and that the remaining terms oc k° yield the results for D2m,2n{t) 
presented in Sec. |Vl 

The next term Z4 requires the introduction of two auxiliary variables ki and K2- Starting from (j22p we have to 
regularize the overlap element 

(ei,-ei,e2,-e2i-6,6,-6,ei> = (ei,-ei,e2,-e2i-e2 + ^^2,e2 + «2,-6 + ^i,Ci + «i) (A2) 

= (27r)4 [5(ei - a - «i)^(-ei + a - «i)^(^2 - 6 - ^2)s{~e2 + 6 - ^2) 

- siCi - a - «i)'5(-ei + a - >^2)sie2 - a - ^^2)<5(-e2 + Ci - '^i) 

- 5(ei - a - 'i2)s{~([ + a - ^1)5(^2 - a - ^i)s{-(2 + a - ^2) 

+ s{e^ - a - 'i2)s{~ei + a - ^2)<5(e^ - a - ^^O'^c-Cs + a - «i) + • • .] • (as) 

Here the dots represent 12 further combinations of 5- functions which lead, in analogy to pop . to terms containing 
integrals restricted to intervals like < a < '*« • These terms in turn yield contributions oc k" with n > 1 which vanish 
when performing the k- integrations. Hence we have not written these terms in (jASp . Now straightforward evaluation 
of the four terms yields 

I / i-oo \ 2 -. „oo 

Zi = -S{-2ki)S{-2k2) U d^\K^iOn --5{-2k^~2k2) j d^ K^{^ + k,)* K^{^ ^ k^)* K^{0^. (A4) 



Appendix B: Calculation of D22 

The first term in the expansion (|2ip which involves kinematical poles is C22{t), which after shifting the rapidities 
in the ket by n reads 

C22{t) = j^ ^ K^iCr K^iO (^', ~e\ a\-C + K,^ + k) e2A^"(coshc'-coshO. (Bi) 

We decompose the form factor into its connected and disconnected pieces using (pS)) 

-2TT 5{S.' - e - «) /(-^' + iTT + ir/, -e + «) 

-2^ <5(-^' + ^ - k) /(e' + iTT + iyy, ^ + «;) 

+2^ <5(C' + ^-k) fi-e + i^ + iT?, e + '^) 

+27r S{-^' - ^ - k) /(e' + ITT + i7?, -e + k) 

+/(C' + ITT + i??i, -e' + iTT + ir/2, -e + K, e + '«) (B2) 

with r],rji — > 0"*". Insertion of (jB2p into (jBip yields three different types of terms, which we denote by C22, C22(t), 
and C|2(t) respectively. 
The first line simply gives 



/•OO 

C2°2 = a(5(-2«;) / d$|i^q(0| 
Jo 



The second term is obtained from the second to fifth lines in (|B2[) . which yield 

C^^it) = ia COth ^^^ /"" ^ Xq(^ + «)* i^q(0 e2Mit[cosh(C+K)-cosh«]_ 

Hereby we have already omitted terms of the form 



(B3) 



(B4) 



ia COth "^^^Jl r g A'q(K - 0* K^iO e2Mit[cosh(K-0-coshC] ^ ^2 ^gg) 



which vanishes in the K-regularization scheme due to 



/ dKP{K) k" -> for L ^ 00, n> 1. 



We can further analyze (IB4I) by expanding the integrand up to 0{k) 

K^ii + «)* K^iO e2Mit[co.h(c+«)-co.hC] ^ |^^(^)|2 ^ ^^^(^^ c«^(^^ ^ 2Mit«;|Xq(e)|' sinh^ + Oi^'). (B6) 

The contributions from the second and third term vanish as they are antisymmetric under ^ — > — ^, thus we arrive at 

C',,{t) = 2ia coth^^J^^ g \K,m'. (B7) 

Finally, the sixth line yields 

Chit) = a^°°|^irq(e')*ifq(C)tanhC'tanhee^*^«(-^h«'— ^''«) 

X clh ^' + ^-- + -^^^ COth ^'-^-- + ^^^ COth ^'-^ + --'^^ COth ^:±i±^l^^. (B8) 

2 2 2 2 ^ ^ 

In order to isolate the singularities we may shift the ^'-contour to the upper half plane or the ^-contour to the lower 
half plane. Doing so we pick up contributions from the poles at ^' — £, — K + iT]2, C — —£, — K + ir]2 and ^ = ^' + k — ir/2, 
^ = — ^' + K — ijyi, respectively, and we obtain 



^2^2(0 = amc r^ f ^ K^iO i^q(0 tanhe' tanh^ coth^ ^^ coth^ ^-tl g2Mit(coshc'-cosh€) (gg) 

Jo 27r J 2-K 2 2 
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2 ./n 27r 



.iacoth^^^ii^i^^ / ^ 7^,(0 e-^*^'--'^? 



X 



{e(e - k) i^q(e - k) e2^-^"™>^h(5-«) ^ Q(_^ _ ^) ^^(__^ _ ^^ g2Mit cosh(5+«) I (g^Q) 



2 ,/n 27r 



X 



{e(^ + n) i^q(e + k) e-2Mitcosh(«+«) _ Q(_^ ^ ^) ^^(_^ ^ ^-) g-2Mitcosh(C-«)| ^ (gj^) 



where the path 7_ Hes in the lower half plane and was explicitly defined in ([59|) . Expanding (jBlOp and (|Blip in k 
yields 

-2lCrCOth ^ / ^l^q(OI ^KCOth ^ / d^|Xq(0| smh^, (B12) 



where we have used 



dK 1 » oo 

^^^<i(0^(0=o^<i(^)\ -0. (B13) 







/o d^ 2 

Now the first term in (jB12p cancels C\2{i) ^nd we arrive at the final result 

C22{t)=aS{-2K) de|ifq(Or / dCl^q(OI'sinhe + i?^2W = '^^2-ari + i?^2(i), (B14) 

Jo T^ Jo 

where Z2 , F and D22 (t) are defined by (lAip , p7|) , and (jB9p or p8)) respectively. We stress that the first term exactly 
cancels the product Z2C00 = aZ2, while the second term equals the second term in the expansion of tre"'"*. 

Appendix C: Calculation of D42 and D24 

We restrict ourselves to Ci2{t) as C2i{t) — Ci2{t)* ■ For the calculation of C/^it) we follow the same steps as above: 
(i) shift the rapidities in the ket by the auxiliary variable |— ^,^) — > |— ^ + k,^ + k), (ii) analytically continue the 
resulting form factor using ([25)) . (iii) evaluate the terms by shifting the contours of integration for ^^ and ^2 to the 
upper half plane, and (iv) expand the result in k up to 0(k°). In particular one can show that all ill-defined terms 
oc 1/k [see JBTj for a similar term in 0{KT)\ cancel each other. After straightforward calculation we obtain 

CA2{t) = S{-2K)C2o{t) d^\K^{0\'-rtC2o{t) + a -^i^qCO' tanh^e^^^'^^^^^ (CI) 

Jo Jo 2-K 

~2ia f ^ r^k^ie)\K^iO\'t^nhe 
Jy+ 2tt Jo 2tt 

X fcoth ^— ^ - tanh i— ^ - coth ^-^1 + tanh ^±^\ g2Mitcoshc' (^2) 

V 2 2 2 2 ) ^ ' 

+ U,,^r^ ^.(^) tanhe n ^.(^0 tanhe: fjcoth^ ^ coth^ ^ 

X tanh^ ^!lzA tanh^ ^^±^2 g2Mit(cosh5;+coshc^-cosh«) ^^3) 

2 2 ^ \ J 

where the path 7+ lies in the upper half plane and is explicitly defined by {0 < 4>o < 7r/4) 

^+^')^i(.-0o)+i0o, cbo<s<oo. (C4) 

The first term in (|C1|) cancels against the product of Z2 C20 (t) , while the second term corresponds to the second term 
in the expansion of C2o(i) e~'"*. All other terms constitute sub- leading corrections to (|i5)). 
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Appendix D: Calculation of D44 

The calculation follows the same steps as outlined above. The main difference is that we have to introduce two 
auxiliary variables, i.e. the form factor in Ci4{t) becomes 

(ei,-e;,C2,-e2kh6,6,-ei,ei)-^(ci,-ci,e2,-^2ki-6 + «2,6 + «2,-a + «i,ei + «i), (di) 

which is analytically continued using (j25p . After a tedious but straightforward evaluation of the resulting terms up 
to 0{k1, K2) one can explicitly show that all ill-defined terms 

oc coth , c>c coth , ex coth (US) 

cancel each other, where we have employed for example 

/ dm dK2 < 4 coth ^""^ ~ '° P{ki) P{k2) oc L^~"~" ^ for m + n>2. (D3) 

The final result is 

C44(t)=fTZ4 + i(<5(-2^i) + <5(-2K2)) (D'^2{t) ~ ^^i) j d^\K^{0\^ + ^{Ttf - T t D'^^{t) (D4) 

Jo ^T^ 

X J ^K.ia? (coth^ ^ - tanh^ ^ - coth^ ^-^ + tanh^ ^) (D5) 

+2a9\t r^ [ ^ KaiO^ KJO tanhf' tanh^ coth^ ^-^ coth^ ^^ g2Mit(coshc'-coshO pg) 
Jo 27r J^_ 27r 2 2 

^^^J %r^J |f^q(0^q(ei)'i^q(6)tanhe'tanh6coth2^^coth2^^ 



'7+ 

X 



tanh tanh coth ^ h coth ■ 



tanh tanh coth h coth 

2 2 2 2 



2 

2Mit(cosh{'~cosh52) (T)7] 



f TTW f f%n^.(^0^.«^)t--hatanh6 n-th^^coth^^ 

X tanh^ ^l^i^ tanh^ ^^±^ tanh^ ^^^ tanh^ ^1±^ g2Afiti:acosh?;-coshC0, (d8) 
2 2 2 2 ' V ; 



where Z4, T and £'22(0 ^^"2 defined in (|A4p . ([57)1 and (IB9p respectively. The paths 7± are defined in ([5^ and (|C4p . 
The leading contributions are given by the first line (ID4p , the second line (jDSp yields the time-independent term Z544 , 
and (jD6p - (jD8p constitute sub-leading corrections that fall off at least as ^ l/{Mt) in the long-time limit. 



Now D44{t) is obtained by [see (|B14I) 

D^iit) = C^iit) - Z2C22it) + [Zi - Zi)a (D9) 

= (s{~2K)-^S{-2K,)-^S{-2K2)]{art^D'22{t)) l^°"d^\K^{0\^ (DIO) 

+ ^(nf-TtD'22{t) + D'44 + ..., (Dll) 

where the dots represent the sub- leading terms (jD6p - (|D8P . Here (jDlOp vanishes due to 

6{-2k) - ^S{-2ki) - ^6{-2k2) = i dK5{-2n) P{k) -]- dKidK2[S{-2Ki) + S{-2k2)]P{ki) P{k2) {B12) 

= f - I / ^^^2 Pi^i2) -^ I dfii P(Ati) = 0, (D13) 

and we arrive at P^ . 
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Appendix E: Leading time dependence of D2m,2m and D2m+2,2m 

The leading behaviour of D2m.2m (t) in the long-time limit is obtained from C2m.2m (t) given in (|2ip by (i) introducing 
the auxiliary variables k^, i = 1,...,to, (ii) regularizing the form factor according to (1251) while keeping only the 
connected piece (i.e. the term with A2 — B2 = $), (iii) evaluating the connected form factor (J14p . and (iv) shifting 
the ^j'-contours to the upper half plane and keeping only the contributions from the poles at ^'^ — ^j — Kj + iO. The 
result after these steps reads 



ncoth^^°°^^i^ni.,(c.-..)*x,te 



a^^ TTcoth ^^i^ r '^^\:':'^^^"' TT KJe^ - tii)* Kj^i) e2Mi*(cosh(C-«.)-coshC.) 



m\ - , ,, , 



fj tanh^''^^+"''"^ tanh ^' ' ^^ ' "' + "^^ coth ^' ' ^^ " "' " "^ + ^° coth ^' ~ ^^^ + "' + "^^ ' ^° . ^^"'^ 



i<j 



Now expanding in the k^'s we obtain 



i°° ^72^ n l^.(^^)l' --he. + • • • = £i (-^ /°° ^^ l^.(Or -nhc) f" + . . . , (E2) 



m! — . 

2 = 1 



where the dots represent terms that grow at most as oc t"^~^ for large times. As the disconnected pieces of C2m.2m{t) 
also grow at most as oc t'^~^ we deduce that the leading time dependence of D2m,2m{t) is given by (|44|. Following 
the same line of argument for C2m+2,2m{t) we arrive at (|^5|) . 



Appendix F: Finite-size regularization 

In our calculation in the infinite-volume system we have regularized the kinematical poles in the form factors using 
a combination of the analytic continuation (|25p together with the K-regularization. An alternative procedure to 
regularize the kinematical poles is to studjsiSii^i^ the system with a finite length L. In this case the Hilbert space 
divides itself into two sectors: the Neveu-Schwarz (NS) sector corresponding to antiperiodic boundary conditions and 
the Ramond (R) sector corresponding to periodic ones. The rapidities in these sectors are quantized according to 

NS : MLsinh^p = 27rp, p G Z + i, (Fl) 

R: MLsmheg = 2Trq, qeZ. (F2) 

The spin operator connects the two sectors, its form factors in the finite system read^ 

m n 

m{Pl,---,Prn\<j\qi,---, qn)ji = S{L)Y[g{^P.) H 9{0qj)F-mn{£.pi,- • ■ , ^p™ l^qi,- ■■,dqj, (F3) 

j=l j=l 

where the function F,„„ is the infinite-volume form factor [see ((14)) and ((27)) ] 



m n 



i^^„(ei,.-.,e™|0i,-.-,^n)=iL("+"^/'J^ n tanh5l_il [] tanh ^^ J] IT '^"th ^1_^ , (F4) 

i,j = l ^-i=l i=lj = l 

the symbol [ J denotes the floor function, and the constant as well as the leg factors are 

S{L) = l + 0{e~^)^ g{e) = ~g{e) = — ^ +0{e-^). (F5) 

V ML cosh 

We stress that due to the quantization of the rapidities the singularities in the form factor ((F4[) are regularized. 

As we consider quenches in the ordered phase which breaks the Z2 invariance the initial state in a system of length 
L has the forroiSiSi 



I*o)l = ^[|*o)ns + I*o)r 



(F6) 
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FIG. 1: D22{t) (black line, left axis) and D|'2(t) + r^t (red line, right axis) for tq = 0.1, L = 50 and iV = 300. For comparison 
we show D22{t) (black stars) and D'22{t) (red stars) as defined in (I36|l and (I38|l respectively. We observe excellent agreement 
between the finite- volume regularization and the infinite- volume results. 



where 



/ 



I*o)ns/r = cxp 



J2 K^{(k)AH-^,)A\(,k 



fceNS/R 

fc>o 



|0) 



NS/R 



(F7) 



/ 



and |0)j,^g /j^ denotes the vacuum state in the corresponding sector. In the regularization introduced by Fioretto and 
Mussardo (see Sec. IVII) the quench matrix Kq{^k) has to be replaced by Krg{£,k), which we will do in the numerical 
evaluations presented in Figs. [T] and [5] The time evolution starting from the initial state (jF6l) now reads 



Ht)) = 



2Ns(^okWl^o)R 
ns(*o|*o)ns+r(*o|*o)r 



— 2-^ '^2m,2nV')- 



m,7i^0 



(F8) 



Straightforward calculation of the norms gives 

oo 
NS/r(*o|*o)nS/R = / . ■^NS/R.2n> 



^NS/R,2 



^NS/R,4 



^ ifq(a)*ifq(6)NS/R(fc,-fc|-',0NS/R- E l^q(^fc)l'' 



A:,ZeNS/R 
fe,/>0 



fcGNS/R 
fc>0 






E i^q(^^)i' -\ Y. i^qfe 



V 



feGNS/R 
fe>0 



/ 



fcGNS/R 
fe>0 



(F9) 
(FIO) 

(Fll) 



We note the analogy to the infinite-volume results (|A1|) and (|A4p respectively. The first non-trivial term in the 
expansion (jFSp reads 



-D22(^) — C'22(i) - - \Zj^^2 + ^R,2) 



CkS) = E E^q(^p)*^q(^'^)Ns(P,-Pkh9,9>Re2*'"(™^'^«''-™'^^«'^ 



pGNS gGR 
p>0 g>0 



^E Y. ^'^^^vT K4l,) 



{ML) 



pGNS gGR 
p>0 g>0 



tanh^p tanh^g 
cosh^p cosh^g 



coth^ 



^?^Ji coth^ ^-^^^ 



^2A'/ii{cosh ^p— cosh^g) 



(F12) 
(F13) 

(F14) 
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FIG. 2: D^i{t) (black line, left axis) and O'^'^it) = D^^it) - {V^tf /2 + r^t(D|'2(t) + V^t) (red line, right axis) for ro = 0.1, 
Z/ = 50 and A^ = 300. For comparison we show Dalt) (black stars) and D44 (red stars) as defined in (|42|) and (|43|) respectively. 
Note that in contrast to -D44 the finite- volume term D44 (i) depends on the time as it also includes corrections which fall off as 
powers of Mt [these corrections are incorporated in the dots in (|42[)] . We observe excellent agreement between the finite- volume 
regularization and the infinite-volume results. 



In the same way a straightforward calculation yields 



^mW 



77^1 E E K^{^v^rK^i^P^rK^{U)K^{U\{ 



A{ML) 



Pi,P2eNS i3i,g2eR 
Pi,P2>0 qi,(;2>0 



tanh^p. tanh^q; 
\^ coshfp. cosh ^5. 



X tanh 
2 



Spi ' t,p 



tanh' 



spi 4p 



tanh 



2 Sgi + S132 



tanh 



2 Sgi sg2 



X 



Jl coth^ klllsL coth^ ^^' „ ^■- e 



2 2 

«£_ 2MitJ2iicosh^p.-cosh^g.) 



i,3 = '^ 



2 ^22(0 (■^NS,2 + ^R,2) + ■^ ■^NS,2 ^R,2 + T (-^NS^ + ^KA) 



Finally, the relaxation rate (HT)) up to ©(i^T^) is given by 



r^ = -- 



2M 

IT 



2M 



^ |Xq(ep)|' tanh^p - ^- E l^qfe)!' tanhe 



pSNS 
p>0 



geR 
g>0 



(F15) 



(F16) 



We have evaluated (jF14|) and (IFISP numerically using the replacement Kq -^ Kt„ for several values of tq, L and N, 
where A'' denotes the UV cut-off for the momentum numbers p and q. The results are shown in Figs. [1] and [51 We 
observe excellent agreement with the results obtained in the infinite volume. 
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